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We have applied the bond valence method to cerium oxides to determine the oxidation states of the 
Ce ion at the various site symmetries of the crystals. The crystals studied include cerium dioxide and 
the two sesquioxides along with some selected intermediate phases which are crystallographically 
well characterized. Our results indicate that cerium dioxide has a mixed-valence ground state 
with an /-electron population on the Ce site of 0.27 while both the A- and C-sesquioxides have 
a nearly pure configuration. The Ce sites in most of the intermediate oxides have non-integral 
valences. Furthermore, many of these valences are different from the values predicted from a naive 
consideration of the stoichiometric valence of the compound. 



I. INTRODUCTION 

Mixed valency of transition metal and rare-earth ions 
in solids and compounds is a question of fundamental 
interest in materials physics, chemisty, and molecular 
biophysics .'i' In 1967, Robin and Day^ published a clas- 
sification scheme for mixed-valence that is still widely 
used today.l^l Class 1 describes systems with two crys- 
tallographic sites that are clearly distinct and and the 
two sites have integral but unequal valence. There is a 
large energy associated with transfer of electrons between 
sites. At the other extreme is Class 3 for which there are 
two sites which are not distinguishable, and one assigns a 
non-integral valence to both sites. The valence electrons 
are delocalised between the two sites. Class 2 is the inter- 
mediate case where the environments of the two sites are 
distinguishable but not very different. The energy associ- 
ated with electron transfer is sufficiently small that it can 
be thermally activated and be associated with significant 
optical absorption in the visible range. On the time scale 
of the vibrations of the atoms the electrons may appear 
to be delocalised. Classes 1 and 2 correspond to what 
Varma^' terms inhomogeneous mixed valence, although, 
perhaps, inhomogeneous integral valence may be more 
appropriate. Class 3 corresponds to homogeneous mixed 
valence. 

In systems which are characterized by homogeneous 
mixed- valence, each ion has the same, noninteger, va- 
lence which is a result of a quantum mechanical super- 
position of two integral valences occuring on each ion 
(See for example Eq. (|2|. Compounds exhibiting this 
type of mixed-valence include, for example, CePda (Ref. 

TmSe (Ref. |6|) and SmBg (Ref. ^ where the va- 
lences of the ions are 3.45, 2.72 and 3.7 for Ce, Tm, 
and Sm respectively. In TmSe, the valence of 2.72 for 
the Tm ion is a result of valence fluctuations of this ion 
between the Tm^''' and Tm'^^ states.^ In contrast, the 
inhomogeneous mixed- valence case involves a mixture of 
different integer valence ions which occupy inequivalent 
lattice sites in a static charge-ordered array. Examples 
of this are provided by Fe304 (Ref. Ej), EU3O4 (Ref. |9|) 
and EU3S4 (Ref. ITD]| . The inverse spinel crystal struc- 



ture of magnetite (Fe304) is considered the classic case 
of inhomogeneous mixed-valence. In this crystal, the 
Fe'^^ ions completely occupy the tetrahedral {A) sublat- 
tice while the octahedral {B) sublattice is equally shared 
between the Fe'^^ and Fe^+ ions and the ionic formula- 



tion is (Fe 



^+) (Fe^+)2 (O^^)^. However, we note that 
this inverse-spinel charge ordering in Fe304 has recently 
been challenged in favour of the normal spinel charge 
structure where the Fe'^"'" ions exclusively occupy all the 
octahedral sites while the Fe^+ ions reside in the tetra- 
hedral sites.l^ The crystal of EU3O4 is a good example 
of the case where ions of different valence strictly occupy 
inequivalent cation sites. It consists of two nonequiva- 
lent Eu sites in which the Eu^+ and Eu'^^ ions occupy 
the eight- and six-coordinated sites respectively giving a 
static charge-ordered array whose ionic formulation may 
be written as (Eu2+) (Eu^+)2 (O^")^ (Ref. ^. 

Oxides of cerium (Ce) appear to exhibit mixed- valence 
characteristics which may help to explain some of their 
properties relevant to their engineering applications. An 
important industrial use of Ce oxides is as anode materi- 
als in high-temperature solid-oxide fuel cells. - For these 
applications, Ce02 (ceria) and Ce203 represent the ex- 
tremal oxidation states in the reversible chemical reac- 
tion, Eq. ([1]): 



Ce02 ^ CeOa-t 



O2 



< y < 0.5 



(1) 



Strong electron-correlation effects may be involved in the 
crystallographic phase transitions which occur between 
Ce02 and Ce203, i.e., the reaction in Eq. ([!]). When 
Ce02 is reduced to the various defective phases, Ce02-y, 
O vacancies are formed in the lattice structure. The 
crystal structure adopted by any such defective phase, 
Ce02-i,, is understood to be the one that provides the 
most favourable energetics for the arrangement of all the 
O vacancies within the structure. In a widely accepted 
view of the microscopic description of O vacancy forma- 
tion and ordering in Ce02-y phases, the two electrons left 
by the O atom when an O vac ancy form s fully localize on 
two of the nearest Ce''+ ions.E^MElil] The locaHzation 
of an electron on a Ce*"*" ion converts it to the slightly 
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larger Ce'^"'" ion with one electron in the 4/ orbital. In 
the reverse process where a defective phase, Ce02-y, is 
oxidized, two 4/ electrons are transferred from the two 
neighbouring Ce^"*" ion sites into the O 2p band. 

This description leads one to expect that the Ce lattice 
sites in the defective Ce02-y phases would consist of a 
mixture of Ce'^"'" and Ce''"'" ions in a static charge-ordered 
array. Thus, useful insight into the microscopic processes 
involved in the reversible chemical reaction, ([T|), could be 
gained from a knowledge of the valences of the Ce ions in 
the defective Ce02~i, phases. However, as we discuss in 
the next section, it turns out that the task of establishing 
the oxidation states of the Ce ions in the lattices of the 
crystal phases involved in ([l]) is very challenging owing 
to valence fluctuations on the Ce ions. 

Here, we report the results of calculations bas ed on a 
simple empirical method, the bond valence modeP^'^^'^, 
to determine the /-electron occupancies in the various Ce 
ion sites for seven of the crystallographic phases involved 
in Eq. ([T]). The crystals studied include j4-Ce203, C- 
Ce203, CesOs, Ce70i2, Ceii02o, CeeOu and Ce02. 

The structure of this paper is as follows. In Section [TTl 
we outline the general problem of mixed-valence in Ce 
oxides from both experimental and theoretical perspec- 
tives. Section Hill describes the bond valence model and 
how we used it to determine the valences of Ce ions in 
the different phases of the Ce oxide crystals. We present 
our results in Section |IV| which we then discuss in light 
of two other models for predicting cationic valences in 
crystals. The main conclusions of this paper are given in 
Section |V| Appendix A lists all the Ce-0 bond lengths 
used in our calculations for Ce-centred polyhedra. 




(»)|/v-> (b)|/y> 

FIG. 1: A schematic of the energy level structure of Ce02 
showing the two important system parameters in the Ander- 
son impurity model: Ae, the energy gap between the two 
states If^p^) and \f^p^) which are mixed by the hybridiza- 
tion, Vpf. For Ce02, the two parameters, Ae and Vpf are com- 
parable and the Ce 4/ level may be close to the O 2p level.'^ 
In (a), the pure /° configuration corresponding to Ce*+ at 
the Ce site is shown. In (b), an electron has hopped from the 
O 2p site to the Ce site where it occupies the 4/ level giving 
the pure configuration, i.e., Ce^+ as shown. An electron 
hole is thus created in the O 2p valence band. The mixing 
between these states due to the hybridization, Vpf, may re- 
sult in mixed-valence at the Ce site. Here, we have set the 
on-site Coulomb repulsion, Uff, (for double occupancy of the 
/ orbital) to infinity. 



II. THE MIXED- VALENCE PROBLEM IN 
CERIUM OXIDES 

In Ce oxides, the number of / electrons on a Ce site, 
Nf, is observed to lie between and slightly above 1.0. 
An isolated Ce atom has the electronic configuration 
[Xe]4/-'5d^6s^. In a crystal of the oxide, the 5d^6s^ elec- 
trons participate in bonding, forming part of the O 2p 
valence band, while the state of the 4/^ electron depends 
on the energy of the associated electronic state within 
the particular crystal structure. If the 4/-^ state lies well 
within the band gap between the empty 5d6s conduction 
band and the O 2p valence band, then the 4/^ electron 
is fully localized on the Ce atom and Nf = 1.0. On the 
other hand, if the 4/^ electron state merges with the va- 
lence band, then the electron is lost to neighbouring O 
atoms so that Nf = thus forming an ionic bond. Inter- 
mediate behaviour with fractional occupancy of the 4/ 
state occurs when the associated state lies close to the va- 
lence band. Occupancies for which Nf > 1.0 are however 
energetically discouraged by the large on-site Coulomb 
repulsion for localized Ce states. 

Fig. [T] illustrates this situation for Ce02 showing the 
relevant energy scales when the system is described by 



the Anderson impurity model.!^ Here, Uff, the on-site 
Coulomb repulsion for the /-orbital on a Ce site is con- 
sidered to be very large. As already noted above, setting 
U f f ^ oo excludes from the ground state wavefunction 
the state |/^p°) which corresponds to Nf = 2.0. Once 
this assumption has been made, then, as shown in Fig. [T| 
the only key parameters of the model become the energy 
gap between the \f^p^) and |/V^) configurations, Ae, 
and the hybridization strength, Vpf. 

Due to the mixing of states shown in Fig. [T| the elec- 
tronic ground state wavefunction of a Ce site in the oxide 
can be written in the general form: 

vI/ = a|/V)+/5|/V) (2) 

where a and fi are constants and 0^ + 0^ = 1, with the 
states assumed to be orthogonal. If for the two states 
|/°p^) and \f^p^), Ae and Vpf are comparable, then the 
nece ssary conditions for valence fluctuation phenomena 
hold Thus one can think of the 4/-level at a given Ce 
site as fluctuating between the /° and /^ configurations. 

Various experimental and theoretical approaches have 
been brought to bear on the problem of mixed valence in 
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Ce oxides with much of the focus concentrated on Ce02 
and to some extent Ce203. There appears to be general 
agreement that the electr onic ground state of Ce203 is 
the configuration.'^^'^^' However, controversy exists on 
the exact nature of the ground state of Ce02 with strong 
arguments in favour of both a pure f'^ configuration and 
a mixed-valence ground state. We briefly review some of 
this interesting debate indicating which approaches have 
consistently arrived at the same conclusions and which 
ones have had mixed interpretations. 

The Bd photoabsorption and photoemission spectra 
analysed by the Anderson impurity model consistently 
reach the c onclusion that C e02 has a mixed valent 
ground state.EaMIMlIIllIlIMl A cluster model has also 
been used to interpret some Sd photoemission spectra 
and it was also concluded that the / st ate is strongly 
mixed- valent in the ground state of Ce02.'^^'^ Although 
the results of 3c? x-ray photoemission spectroscopy (XPS) 
interpreted in the Anderson impurity model have consis- 
tently predicted a mixed- valent ground state for Ce02, 
doubts have been raised about the reliability of assign- 
ing initial state configurations from Ce 3d XPS spectra 
because of the possibility of reduction of the oxide on ex- 
posure to the x-ray radiation.!^ Perhaps, an even more 
difficult challenge in the interpretation of XPS spectra is 
how to trace the ground state configuration of a sample 
from the spectral signatures of the final state of XPS. 
The conventional interpretation is that in Ce203, the 
two-peak structure observed in XPS spectra is a result 
of final state effects whereas Ce02 shows a three-peak 
structure in its XPS spectru m which has been attributed 
to final state eSects0^^ 

X-ray absorption near-edge spectroscopy (XANES) 
has provided a second source of evidence for the mixed 
valence of the ground state of Ce02. The charac- 
teristic double peak structure observed in all XANES 
spect ra is considered a signature for mixed valence in 

Q(^Q^ |30I31I32I33I34I35I36I37I 

In contrast to the 3d XPS and XANES data discussed 
above, the 4d-4/ photoabsorption spectra of Ce02 sup- 
port the ground state configuration based on their 
close resemblance to La trihalide spectra which have no 
/ electron while they differ considerably fr om th e Ce tri- 
halide spectra which have one / electron.EHISa Also m 
favour of the f'~' ground state configuration are results 
from refiectance spectroscopy in the Ad-Af absorption 
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region. 

Among the methods which have produced an am- 
biguous picture of the ground state of Ce02 are in- 
cluded resonant inelastic x-ray spectroscopy (RIXS), 
bremsstrahlung isochromat spectroscopy (BIS), valence 
band XPS and energy band structure calculations. Ap- 
plication of resonant photoemission to the problem led to 
the conclusion that the grou nd state of Ce02 is mixed- 
valent by some author j^^ ^ while a pure f° configura- 
tion was claimed by others. '^^'^^'^ The spectrum inves- 
tigated by the bremstrahlung isochromat spectroscopy 
(BIS) combined with XPS study shows empty localized 



4/ states in the band gap provid ing evidence in support of 
the pure /" configuration This conclusion was also 
supported by valence band XPS studies .'^^'^ However, 
it was later shown that the valence-band photoemission, 
BIS, and 4d photoabsorption spectra can all be explained 
consistently with ot her core -level spectra by the Ander- 
son impurity model.'^SEEll xhe results of this analysis 
led to the conclusion that Ce02 is mixed- valent in the 
ground state. 

An energy band calculation by the linear augmented 
plane wave (LAPW) method showed a 4/ electron count 
of ~ 0.5 with considerable covalent character.''^ Simi- 
lar results were obtained from a linear muffin-tin orbital 
(LMTO) band structurtP and in a molecular orbital cal- 
culation of a CeOg cluster.'^ On the other hand, band 
structures of Ce02 obtained from density functional the- 
ory (DFT) calculations, both in the local density approx- 
imation (LDA) and generalized gradient approximation 
(GGA), have empty 4/ states above the vale nce band 
supporting the /° ground state configuration .1^3^11111111 
This electronic configuration for Ce02 is widely accepted 
in DFT work largely because it predicts the structural 
properties of the Ce02 crystal with reasonable accuracy. 
However, LDA and GGA do not give the experimen- 
tally observed band structure and structural properties of 
Ce203. Some pragmatic strategies have been adopted to 
augment LDA and GGA so that predictions closer to the 
experimental results could be obtained. Skorodumova et 
al)^ artificially localized the 4/ states in what was called 
the 'core state model' while others have used the LDA 
(GGA) -I- U formalism, where an onsite Coulomb repul- 
sion is incorporated to account for the repulsive energy 
arising from the double occupancy of an / orbital on a 
Ce site.lSa 

The problem of the 4/ states in the ground state con- 
figurations of Ce02 and Ce203 resembles the difficulties 
also encountered in studying these same states in Ce 
metal. Ce metal undergoes an isostructural a-7 phase 
transition with a 15% change in volume.^^ The average 
/ electron populations observed experimentally are 0.81 
and 0.97 for a- and 7-Ce respectively.'^S Thus, the pure 
Ce metal shows signiflcant mixed valence in the 7 phase. 

III. THE BOND VALENCE MODEL 

The bond valence model (BVM), which is a general- 
ization of Pauling's original concept of the electrostatic 
valence principle , has been reviewed extensively in the 
literaturelSIiZIiSl A quantum chemical justification of the 
model has been discussed by Mohri.'^The model defines 
the relationship between bond valences, s, and the cor- 
responding atomic valences, V, through the equations: 

= E^.j (3) 
j 

where i and j refer to different atoms. This relationship 
is called the valence-sum rule. 
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There is a well-defined empirical relationship between 
bond valences as defined in ^ and bond lengths in a 
given coordination polyhedron and it is this functional 
relationship which makes the BVM quantitatively useful. 
The relationship is monotonic, and over the small range 
in which most bonds are found, it can be approximated 
byP 



exp 



^0 — Ri, 
B 



(4) 



Here Rq and B are fitted parameters, Rq being the bond 
length of a bond of unit valence. For a specific bond, 
these parameters are determined by fitting the measured 
lengths of bonds in a wide range of compounds by en- 
forcing the valence-sum rule (|3|. It has been shown that 
for most bonds, B can actually be set to be 0.37 A which 
reduces the bond valence model Q to a one-parameter 
model.^^ In that case, for each structure where a central 
atom is bonded only to atoms of a particular species, the 
Rq parameter is then obtained by combining the valence- 
sum rule ^ and Q, i.e: 



Rq 



3 3 



B 



exp I — I exp 



Which can be rewritten in the form: 



B 



(5) 



(6) 



The BVM is applicable to bipartite crystals which, in our 
case, means that only Ce-0 bonds can exist in any given 
crystal. In all cases, the bipartite requirement is satisfied 
for Ce polyhedra considered below. However this is not 
true of all O polyhedra as in a few cases, e.g., in CegOn, 
O atoms are included in the coordination polyhedron of 
an O atom. Nevertheless, this does not concern us as all 
our calculations are performed on cation-centred polyhe- 
dra. A second issue arises from how to precisely define 
the coordination number for some of the Ce sites. We 
adopt Brown's definition of the coordination number as 
the number of atoms (O, in this case) to which a cen- 
tral atom (Ce, in this case) is bonded.'^ The operational 
meaning of this definition was given by Altermatt and 
Brown;^" a bond exists between a cation and an anion if 
its experimental bond valence is larger than 0.04 x the 
cation valence. 

Various authors have determined values for the param- 
eters Rq and B in Eq. Q for Ce-0 bonds from an 
analysis of measured bon d len gths in both organi(pii£2l 
and inorganic compounds .'^^'^ It has been reported that 
Ce-0 bonds in inorganic compounds are longer than in 
metal-organic coordination compounds and so the bond 



parameters are larger!^. Although Zocchi has derived 
detailed parameters for the Ce-0 bonds in iorganic com- 
pounds explicitly giving the dependence on coordination 
number^, his so-called method of intercepts'*'^ does not 
provide a clear physical basis for the calculation. In view 
of these considerations, we selected the parameterization 
of the BVM by Brese and O'Keeffe.ES! 

In this parameterization, B = 0.37 A with Rq — 
2.151 A and Rq = 2.028 A for Ce^+ and Ce^+ ions re- 
spectively. This means that, by Brown's criterion for 
the coordination sphere, only O atoms within the critical 
distances Rc = 2.905 A and Rc = 2.746 A are included 
within the coordination polyhedron for Ce^+ and Ce''"'" 
sites respectively. 

To perform bond valence calculations, the only input 
required is the bond length data of the Ce-0 bonds in 
the respective coordination polyhedra of the various ox- 
ides. For this, we obtained crystallographic data from the 
sources listed in Table |Tj For each oxide, all the distinct 
coordination spheres for both Ce and O atoms were iden- 
tified. Fig. [2] illustrates the procedure using the example 
of Ce70i2 where we have shown only one of the two O 
atom polyhedra in the crystal of this oxide. The Ce- 
O bond distances in each polyhedron were then obtained 
from the crystallographic data using CrystalMaker.lSSlWe 
have listed the Ce-0 bond length data in Appendix A. 
We note here two oxides, C-Ce203 and CegOn, for which 
we were not able to obtain full crystallographic data. C- 
Ce203 has not been observed experimentally and the lat- 
tice constant used here, 11.22 A, was derived by Eyring® 
and also by Tsunekawa et al. We could not find posi- 
tional parameters for the crystal of CegOn and we used 
the positional parameters for PrgOn with which it is 
isostructural.'^ 



IV. RESULTS AND DISCUSSION 

We have noticed from the results of Roulhac!^ that the 
parameter Rq varies approximately linearly with the ox- 
idation state of the Ce atom between Rq — 2.121(13) A 
(Ce3+) and Rq = 2.068(12) A (Ce^"^). We exploit this re- 
lationship between Rq and the Ce oxidation state to cal- 
culate the / occupancies of the Ce sites self-consistently 
in the following way. As already mentioned, the valence 
at each Ce site fluctuates between the /° (Ce*+) and 
(Ce^"*") configurations. Let x be the probability that a Ce 
ion is in the configuration (0 < x < 1), then, the corre- 
sponding probability of the configuration at the same 
site is given by 1 — x. Thus, for such a mixed- valence Ce 
site, the correponding value of the parameter Rq is then 
given by linear interpolation between these values: 



Rq = R^x + i?^ (1 - x) 



(7) 



As already mentioned, we have used i?!] = 2.151 A and 
Rl = 2.028 A, i.e., Rq parameters for Ce^+ and Ce'*"'" 
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FIG. 2: Identification of tlie coordination polyfiedra in 
Ce70i2 for bond valence calculations. Atom labels refer to 
the distinct sites in the crystal, (a) The unit cell of Ce70i2 in 
the hexagonal crystal lattice. The space group is R3 and the 
unit cell consists of three formula units, (b) The polyhedron 
of the Ce(2) site which is seven-coordinated and of triclinic 
symmetry and (c) that of the the six-coordinated Ce(l) site 
of Ss symmetry. There are two distinct polyhedra for the O 
atoms and one of these, the 0(1) site is shown in (d). In this 
polyhedron, the O atom has a coordination of four in triclinic 
symmetry. The colours of the atoms are: Ce - black and O - 
violet. Images generated in CrystalMakeM 



states respectively.'^ From (|7| and Q , we have: 



f [R'ix + R'iil~x)]-R., \ 
Si] [x) = exp I I (8) 



It then follows that, for a Ce site z, Vi {x) = 4 — x = 
^Sij (a;), which then with Nf — x, leads to the required 
j 

self-consistent equations: 



Eexp 



[R'^x + i?2 (1 - x)] - R, 
B 



(9) 



We solved ^ for x = Nf to get the / occupancies of 
the various (Je sites in the seven crystals studied. The 
results of applying Eq. ^ to the seven oxides of cerium 
selected for this study are given in Table |l] 

For Ce02, all O atoms are symmetry equivalent (as are 
all Ce atoms). For this simple case, we can write down an 
equation describing the polyhedron centred on given O, 
equivalent to ([9|, and use it to determine the valence of 
the O atom self-consistently. The valence on the oxygen 



atom of Ce02 is given by: 

4 4 

Vi (x) = -1.5a; - 2(1 - x) = J^^jt (x) = -J^^ij (x) 



1=1 1=1 



4 [R'ix + R'i{l~x)]~R,, , 
- E exp ) (10) 



4 [ROx + Ri{l~x)]~R,, 
a; = 4 — 2 2^ exp 



i=l 



B 



(11) 



Evaluation of ( 11 ) gives x — 0.268 from where it follows 



that the O valence is —1.87 which is consistent with the 
Ce valence already calculated above. 

In the general case, in which there are inequivalent O 
atoms, we cannot determine the valence of a given O site 
in this way. That is, while the valence of each Ce atom 
can be determined by considering the polyhedron consist- 
ing of this central Ce and its nearest neighbors, the same 
is not generally true of an O atom. This is because we 
have assumed that the parameter Rq, a property of a Ce- 
O bond, depends on the valence of the Ce atom but not 
on the valence of the O atom. This assumption is consis- 
tent with the observation that the structures of different 
Ce-oxides with the same Ce-valence can be adequately 
described by (|4| using a single value for Rq (Ref.79|. 

Locock and Burns have defined a measure of the vari- 
ation in bond lengths of the bonds included in a coordi- 
nation polyhedron as follows JsHl 

var^aUon{%) = l^'^^^-^a. - 6ond^.„| 
bondavg 

where the subscripts max, min, and avg refer to the max- 
imum, minimum and average bond lengths in the poly- 
hedron. We have calculated the % variation for all the 
polyhedra in this study and the results are given in Table 
III 

We now give a few remarks about the accuracy of the 
bond valence method. Brown has estimated that bond 
valence sums have an accuracy of 0.05 v.u.^^ The main 
error in the method comes from the fitted parameters 
of the model {B, Rq). It has been noted that bond va- 
lence parameters which overestimate valences of strong 
bonds while underestimating those weak bonds will give 
bond valence sums which are too high in low coordination 
polyhedra and too low in the case of high coordination.lSll 
We have estimated that with the following uncertain- 
ties in the model parameters and bond length data: Rq 
(±0.01 A), R,j (±0.01 A)58 and B (0.037A)i^ the uncer- 
tainties in the bond valence sums for Ce02 and A— Ce203 
are ±0.13 vu and ±0.12 v.u respectively which is in agree- 
ment with Ref. ^82, 

In Fig. |3] we have plotted the Ce site / occupancies, 
Nf, calculated from Eq. ^ . The results are plotted 
for increasing average degree of oxidation of the Ce ion. 
For each crystal, the results are reported according to 
the site symmetries of the Ce sites which range from the 
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TABLE I: Site valencies of Ce at different sites in cerium oxides (Ce02-t/) calculated from our self-consistent bond valence 
method. 



Oxide 


y 


Ref. 


Ce site 


Coord. 


No. Site Sym- 
metry 


No. of Sites Site Bond Nf 
in Unit Cell Valence 

Sum (v.u) 


% Varia- 
tion 




0.50 




Ce 


7 


Co 


2 


2.97 


1.03 


14 




0.50 


170171 1791 


Ce(l) 


5 




8 


2.98 


1.02 


0.0 








Ce(2) 


6 




24 


2.98 


1.02 


4 




0.33 


I70I72I7'5I 


Ce(l) 


8 




8 


3.33 


0.67 


8.5 








Ce(2) 


8 


^2 


24 


3.37 


0.63 


8.7 




0.29 


1741 751 


Ce(l) 


5 


^6 


3 


3.67 


0.33 


0.0 








Ce(2) 


7 




18 


3.21 


0.79 


7.0 


v_^cil W20 


0.18 


I70|7A| 


Ce(l) 


8 




2 


3.08 


0.92 


13 








Ce(2) 


8 




2 


3.06 


0.94 


7.6 








Ce(3) 


7 


J 


2 


3.46 


0.54 


7.8 








Ce(4) 


7 




2 


3.58 


0.42 


12.5 








Ce(5) 


7 




2 


3.68 


0.32 


8.2 








Ce(6) 


7 


- 


2 


3.76 


0.24 


9.4 


CesOii 


0.17 


I(i8l71l77l 


Ce(l) 


6 




4 


3.22 


0.78 


30.9 








Ce(2) 


7 




4 


3.75 


0.25 


26.7 








Ce(3) 


7 




4 


3.10 


0.90 


15.1 








Ce(4) 


5 




4 


3.62 


0.38 


37.4 








Ce(5) 


6 




4 


2.93 


1.07 


24.0 








Ce(6) 


6 




4 


3.22 


0.78 


30.9 


CeOa 


0.00 


I68I71I78I 


Ce 


8 


Oh 


4 


3.73 


0.27 


0.0 



lowest triclinic sites («) to the most symmetrical octahe- 
dral sites (Ofi)- The dashed line labelled Ce203 refers to 
both the A- and C-phases of this oxide. For comparison, 
we have included the solid straight line which represents 
Nf values calculated from the formula units by requiring 
charge balance of the formula unit, a valence of —2 for all 
O ions and an even distribution of the valence among all 
the Ce ions in a formula unit. We will call this method of 
calculating Nf the homogeneous mixed-valence approx- 
imation (HMA) as it assigns the same valence to each 
cationic site regardless of the specific site properties. We 
contrast this method with another simple approach for 
predicting cationic valences which we will call the inho- 
mogeneous mixed- valence approximation (IHMA). The 
IHMA is based on the requirement that all valences must 
be integral but their exact assignment in the crystal lat- 
tice does not necessarily have to satisfy constraints which 
may arise from specific differences in local site symmetry. 
Thus, this approximation considers the crystal lattice to 
be a static charge ordered array. 

The results plotted in Fig. |3]show that mixed- valence 
in Ce oxides does not fit nicely into either of the tra- 
ditional classes of mixed- valence, i.e., homogeneous or 
inhomogeneous mixed-valence as originally defined by 
Varma.l^ For the Ce oxides, only Ce02 and Ce203 are 
strictly homogeneous mixed-valent oxides since all the 
Ce sites are symmetrically equivalent and have the same 
oxidation state. The rest of the oxides do not have sym- 




1.6 1.8 
0/Ce raoinr ratio 



FIG. 3: Ce /-level site occupancies in different crystallo- 
graphic phases of the oxides calculated self-consistently in the 
bond valence model. The results are given according to the 
exact point group symmetries of the respective Ce sites in 
each crystallographic phase. As can be seen from the legend, 
the point symmetries of the Ce sites vary from as low as tri- 
clinic (i) to as high as octahedral (Oh), the point group of the 
cube. The straight line represents Nf values obtained from a 
simple Ce valence calculation based on the electroneutrality 
of the formula unit and the assumption of an oxidation state 
of —2 for the O ions in all the oxides. 
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metrically equivalent Ce sites and the oxidation states 
of the individual sites are generally different within the 
error limits of the method. The mixed-valence in the 
crystal is not a result of the averaging of the oxidation 
states at the different sites but a property of the individ- 
ual sites, which means that the mixed-valence is not of 
the simple inhomogeneous type. 

It is interesting and important to compare our results 
to those obtained by other methods. In Table [ll] we have 
reported the results of the Ce site 4/ occupancy, Nf , from 
the literature. We have not been able to obtain results 
for all other oxides except for Ce02 and Ce203. The data 
in Table [ll] shows that while there is good agreement on 
the / occupancy of the Ce site in A-Ce203, there is con- 
siderable variation in the reported results for Ce02 . The 
Nf values for Cc02 vary from 0.00 for LDA and GGA 
calculations to 0.60 for Ce 3d core XPS. We now discuss 
our results in detail for the individual oxides. The format 
of our discussion is as follows: for each oxide, we compare 
the predicted valences from the BVM to those obtainded 
from HMA and IHMA. We then indicate whether the 
particular oxide exhibits mixed- valence or not. We sum 
up this section by relating our results to the conventional 
model of vacancy formation and ordering in bulk Ce ox- 
ides and the possible role of electron correlation in these 
oxides. 



A. ^-CeaOg 

The bond valence method gives Nf = 1.0 for this ox- 
ide. As can be seen from Table |ll] the bond valence 
method gives iV^ values which accord well with results 
from the other methods for this oxide. In addition, both 
the HMA (Fig. |3| and IHMA also predict iVj = 1.0 and 
so all the methods are in agreement. Compared to Ce02, 
Table |ll] shows that the variation in bond lengths in the 
Ce polyhedron of A-Ce203 is relatively high at about 
14%. However, the result of Nf obtained in this calcu- 
lation which is in good agreement with other methods 
may suggest that bond lengths distortions of this mag- 
nitude may have no significant role in the bond valence 
model. We conclude that mixed- valence in A-Ce203 is 
neghgible. 



B. C-CeaOa 

The results are similar to those for A-Ce203 just de- 
scribed above in both the Sq and C2 site symmetries of 
C-Ce203. However, structurally, the two sesquioxides 
are very different. The ^-sesquioxide has a hexagonal 
Bravais lattice (space group P3ml) with a lattice con- 
stants a = 3.891 A and b = 6.059 A and only one Ce site 
of symmetry. In contrast, the C-sesquioxide has a 
cubic fluorite structure with two O vacancies along the 
body and face diagonals (space group /a3) giving the 
Sq and C2 site symmetries for the Ce ions respectively. 



That the Ce sites in A-Ce203 and C-Ce203 still have the 
same Nf value despite the differences in point symmetry 
appears to suggest that site symmetry may have no sig- 
nificant influence on the valence of a given site in Ce203. 
Compared to the other oxides included in this study, the 
only other Ce site with a comparable Nf value is the 
triclinic Ce(5) site in CcgOn. 



C. CeaOs 

The results in Table |l] indicate that both Ce sites 
in this crystal have comparable valences with Nf val- 
ues of 0.67 and 0.63 for the Se and C2 sites respec- 
tively. Again, this result highlights the earlier obser- 
vation that site symmetry may have no significant role 
in site valences as the two Nf values of these two sites 
are comparable. The HMA predicts that Nf = 0.67 
for this crystal which compares favourably with the 
BVM while the IHMA requires the ionic formulation 
[{Ce^+)^Ce^+] (O^-)^ for Ce305. Thus, the IHMA in- 
dicates that in a static charge-ordered array of integral 
valences, the Ce305 lattice has twice as many Ce'^"'" ions 
as there are Ce^"*" ions. Although the Ce(l) and Ce(2) 
sites are of different point group symmetries, they are 
of comparable polyhedron sizes with average Ce-0 dis- 
tances of 2.444 A and 2.435 A respectively. When one 
considers how to distribute the Ce'^"'" and Ce"*"*" ions be- 
tween these sites, there is a conceptual difficulty. Firstly, 
in typical ionic crystals of Ce, the average radii of the ions 
in an 8-coordinate environment are 1.28 A and 1.11 A for 
Ce^+ and Ce^'^ respectively.'^ This 13% difference in the 
crystal radii of the ions is significantly different from the 
0.4% difference in the average sizes of the coordination 
polyhedra. Secondly, the ratio of the Ce(l) and Ce(2) 
sites does not match that of the ions as given in the ionic 
formulation above. Thus, if one were to argue that the 
smaller Ce^"*" ion will show a slight preference for the 
smaller Ce(2) site, then all the Ce*"*" ions would occupy 
the Ce(2) sites. However, since there are almost twice as 
many Ce(2) sites as there are Ce^+ ions, the remainder 
of these sites will be occupied by the Ce'^"'" ions and the 
latter, will, in addition, occupy all the Ce(l) sites. This, 
of course, results in a situation where Ce^"*" and Ce'*+ 
occupy the same type of site in the crystal. This logical 
inconsistency does not arise when the BVM is applied 
and we conclude that Ce305 is a homogeneous mixed- 
valence compound. 



D. CerOia 

Both the Ce(l) and Ce(2) sites are predicted in the 
BVM to be mixed-valent with the more symmetrical 
Ce(l) site closer to Ce^"*" (3.67 v.u) and the triclinic 
Ce(2) site closer to Ce^+ (3.21 v.u) whereas the HMA 
gives a valence of 3.4 v.u. The ionic formulation for the 
IHMA is [(Ce^+)^ (Ce*+)3] (O"^)^^ 
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TABLE II: Nf values for the Ce ion in cerium oxides determined by different methods 



Compound 


Nf 


Method 


Remarks 


Ref. 


CeOa 


0.3 ±0.1 


Bond valence 


From crystal structure data 


This work 




0.60 


Ce 3d core XPS 


Bulk measurement 


[27] 




0.45 


XANES 


Bulk measurement 


[H3 




0.50 


Ce 3d core XPS 


Bulk measurement 


[I9] 




0.05 


Optical reflectance 


Bulk, measurement 


HI 




0.2 - 0.4 


LDA +U 


Bulk, calculation 






0.50 


LAPW Xc, 


Bulk, calculation 


m 




0.40 


Ce 3d XAS 


Bulk, calculation 


[86] 




0.20 


HSE 


Bulk, calculation 






0.00 


LDA and GGA 


Bulk, calculation 


M 


CeaOg 


1.0 ±0.1 


Bond valence 


From crystal structure data 


This work 




1.0 


Ce 3d core XPS 


Bulk measurement 






1.0 


LDA and GGA^° 


Bulk, calculation. 


M 




1.0 


HSE 


Bulk, calculation 


M 


Ce2Zr207.5 


0.5 


EELS (Ce-M4,5) 


Bulk measurement, ceria- 


m 



zirconia solid solution 



^''Both the LDA and GGA calculations were performed by artificially localizing the 4/ states, the so called core state model. Hence, 
A'^ = 1 was actually assumed in the calculation. 



[{Ce^+)^^Ce^+] (0-2)3g for CeyOia. The Ce(l) site is 
smaller than the Ce(2) site and it has been suggested 
that since the lower coordination at the Ce(l) site would 
require a smaller cation, then Ce^"*" ions are expected to 
occupy all three of these sites. The remaining six Ce'*+ 
and the twelve Ce^"'' then occupy the Ce(2) which leads 
to a similar problem as already noted for CeaOs above. 
When the BVM and IHMA results are compared, we 
notice that the BVM predicts that there are two distinct 
valences for the Ce ions, one closer to Ce'^"'" and the other 
closer to Ce'*+but without the inconsistency of assigning 
different valences to the same type of Ce site. These 
results indicate mixed- valence in Ce70i2. 



E. CeiiOao 

This oxide has the lowest symmetry of all the cerium 
oxides included in this study. Based on the the BVM re- 
sults in Table [T] the following approximate assignments 
of site oxidation states can be made: Ce(l) and Ce(2) 
- Ce'^+j Ce(3) and Ce(4) - strongly mixed- valent, Ce(5) 
and Ce(6) - closer to Ce''+. Fig. jsjshows that the HMA 
predicts a valence of 3.6 v.u. which is also mixed- valent. 
The ionic formulation [{Ce^+) ^ {Ce'^+) ^] (O-^)^^ results 
from applying the IHMA to this crystal. Both the BVM 
and the IHMA predict the presence of Ce'^+ sites in this 
crystal and a consistent distribution of the sites is ob- 
tained between the methods since it is expected that, in 
the IHMA ionic formulation, the larger Ce^+ ion will oc- 
cupy the larger eight-coordinated polyhedra Ce(l) and 
Ce(2) which is the same result obtained from the BVM 
analysis. However, the two methods disagree on the as- 



signment of valencies to the remainder of the sites with 
the IHMA assigning them all to Ce'^+ while the BVM 
clearly indicates strong mixed-valence for the Ce(3) and 
Ce(4) sites and less but still significant mixed valence for 
the Ce(5) and Ce(6) sites. The HMA does not appear to 
be a viable proposition for this crystal lattice because of 
its very low symmetry. Again, we find that Ceii02o is a 
mixed-valence compound but one that does not fit into 
either of the traditional classes. 



F. CesOii 

The results in Table |T] obtained from the BVM show 
that all Ce sites in CegOn are mixed-valent with de- 
viations from the nearest integral valences increasing 
from 0.1 v.u. to 0.38 v.u. in the order Ce(5), Ce(3), 
Ce(l), Ce(6), Ce(2) and Ce(4). As shown in Fig [s] 
the HMA predicts mixed-valence for this compound giv- 
ing a valence of 3.7 v.u for the Ce ions. Considered in 
the IHMA, the ionic composition of CegOn should be 
[(Ce''^+)2 (Ce''+)J (O^-)^^ which implies that there are 
twice as many Ce^+ ion sites as there are Ce^"*" sites in a 
unit cell of CcgOn. We notice that the predictions of the 
three methods are very different for this crystal. Table |I] 
shows that this crystal has the most distorted polyhedra 
with % variation in the bond lengths ranging 15% — 37% 
and it is expected that the BVM would reflect this aspect 
of the local site geometries. With these complex local ge- 
ometries, it is not expected that the HMA would give a 
good approximation to the site valences. Comparing the 
IHMA and the BVM, we notice that even if one tenta- 
tively considered both Ce(2) and Ce(4) to be Ce^"'" sites 



9 



and the rest Ce^+, then one gets twice as many Ce^+ 
sites as there are Ce^"*" sites. This charge-ordering is the 
reverse of that predicted by the IHMA ionic formulation 
which clearly highlights the disagreement between the 
methods. 



G. CeOa 

Our value of Nf = 0.27 for Ce02 is comparable to the 
result obtained by Castleton et al^ from an LDA + U 
calculation as shown in Table |TT1 We can see from Ta- 
ble [n] that results from 3d core-level spectroscopy tend 
to give relatively high Nf vaues {Nf > 0.40). In Ce02, 
the Ce site is in a symmetric polyhedron with all the 
Ce-0 bonds equal in length and therefore the bond va- 
lence model is expected to perform well for this oxide. 
As can be seen from Fig. [3j the HMA predicts that the 
Ce ion is in a pure /° configuration in Cc02. The same 
result is obtained from the IHMA and, thus, both meth- 
ods contradict the BVM. From our result, wc conclude 
that Ce02 is a mixed-valent compound. 

It would be informative to compare the bond valence 
sum results obtained here to those for their Pr oxide 
counterparts. However, we have not found any detailed 
study of the bond valence sums of Pr oxides for a mean- 
ingful comparison to be made here. We only found a 
bond valence calculation performed for the Pr-0 bond in 
the high-Tc superconductor PrYBa2Cu307 where the Pr 
was mixed-valent with a valence of 3.4 v.u.l^ 



V. CONCLUSION 

Since we did not obtain integral Nf values for most of 
the Ce oxide sites, our results suggest that mixed valence 



is an essential feature of Ce in its pure oxide phases. We 
have also shown that for a given crystallographic phase, 
several different oxidation states may exist for Ce sites 
which belong to the same point group. Thus, it appears 
that valence fluctuations depend on the exact coordina- 
tion geometry of a Ce site which suggests that valence 
fluctuation is a 'local' property of the Ce site. The point 
group symmetry of a Ce site does not appear to play a 
signiflcant role in determining the valence state of that 
site. We have also noted that the bond valence method 
eliminates some of the conceptual difficulties which arise 
when one attempts to distribute Ce ions of different inte- 
gral valences in sites of the same type in a given crystal. 
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APPENDIX A: BOND LENGTHS IN THE 

VARIOUS CERIUM COORDINATION 
POLYHEDRA OF THE CERIUM OXIDES 

The calculated Ce-0 bond distances in the various Ce 
coordination polyhedra of the Ce oxide crystals are listed 
in Table HTH 
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TABLE III: Bond length data used in the bond valence calculations of the Ce oxides. Crystallographic data was obtained from 
the references as cited in the table and the Ce-0 bond lengths were calculated from that data in CrystalMaker® !^ The O 
atoms, 0(i) (i = 1 — 8) in each polyhedron may be equivalent or distinct O sites in a given crystal. All bond lengths are in A. 
Estimated errors in bond lengths are approximately 0.01 A. 



Oxido 


Rcf. 


Ce Site 






Bond Lenf 


;ths in Ce Coordination Polyhedron 












0(1) 


0(2) 


0(3) 


0(4) 


0(5) 


0(6) 


0(7) 


0(8) 


A-CeaOa 


I69IV0I 


Ce 


2.339 


2.339 


2.339 


2.434 


2.694 


2.694 


2.694 




C-CeaOs 


I70I71I72I 


Ce(l) 


2.416 


2.416 


2.416 


2.416 


2.416 


2.416 










Ce(2) 


2.375 


2.375 


2.403 


2.403 


2.469 


2.469 






CeaOs 


I70I72I73I 


Ce(l) 


2.393 


2.393 


2.393 


2.393 


2.393 


2.393 


2.600 


2.600 






Ce(2) 


2.352 


2.352 


2.380 


2.380 


2.445 


2.445 


2.563 


2.563 


Ce70i2 


|74|75| 


Ce(l) 


2.249 


2.249 


2.249 


2.249 


2.249 


2.249 










Ce(2) 


2.405 


2.392 


2.652 


2.307 


2.367 


2.380 


2.478 




CeiiOao 


|7:^|76| 


Ce(l) 


2.489 


2.474 


2.379 


2.489 


2.474 


2.379 


2.701 


2.701 






Ce(2) 


2.423 


2.407 


2.446 


2.624 


2.589 


2.444 


2.539 


2.597 






Ce(3) 


2.311 


2.405 


2.263 


2.333 


2.447 


2.385 


2.379 








Ce(4) 


2.251 


2.309 


2.361 


2.544 


2.267 


2.279 


2.356 








Ce(5) 


2.429 


2.254 


2.268 


2.342 


2.240 


2.358 


2.281 








Ce(6) 


2.214 


2.280 


2.227 


2.315 


2.306 


2.274 


2.429 




CeeOii 


I68I71I77I 


Ce(l) 


2.180 


2.241 


2.245 


2.490 


2.560 


2.618 










Ce(2) 


2.189 


2.230 


2.236 


2.252 


2.516 


2.669 


2.834 








Ce(3) 


2.260 


2.324 


2.377 


2.441 


2.593 


2.603 


2.698 








Ce(4) 


2.064 


2.143 


2.147 


2.189 


2.670 












Ce(5) 


2.235 


2.285 


2.347 


2.482 


2.723 


2.830 










Ce(6) 


2.180 


2.241 


2.245 


2.490 


2.590 


2.618 






CeOa 


I68I71I78I 


Ce 


2.434 


2.434 


2.434 


2.434 


2.434 


2.434 


2.434 


2.434 



